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I. Introduction 


Observations in the soft part of the spectrum indicate 

that most flares produce a hot quasi-thermal plasma at coronal heights 
with maximum temperatures of order 3 X 10 K. The dominant cooling 
processes are believed to fee radlstibn and conduction to the chromosphere 
(Culhane et al., 1970). Simple estimates based on a specific flare 
model (Sturrock, I968) indicate that, in this temperature regime, 
conduction cooling is more important than radiation. This conclusion 
is supported by detailed calculations by Strauss and Papagiannls (1971- ) 
on the same flare model, and by data analysis carried out by Moore and 
Datlowe ( 1975 ) on 17 small flares. 

- If conduction is Indeed the dominant cooling mechanism, then 
oagnstic-field geometry will have a strong influence on the cooling 
rate. Optical observations (see, for instance, Fisher, I 97 I) show 
that the coronal post-flare plasma typically has a loop-like structure, 
presumably due to the influence of magnetic fields. Rust: aiui Bar (1973) 
have deduced from their observations that t^e field in post-flare 
regions is well represented by a potential fields, especially at large 
heights. In the case li^ilch they investigated (the flare of August 7> 
1972), the field strength at the top of the loops was of order 50 gauss 
while that at the base of the loops was over lOOO gauss. This shows 
that there is a large change in cross section of a flux tube along its 
length, a factor which is not taken into account in calculations mentioned 
above. 

It therefore seems desirable td study the^evolutlon of a h 
post-flare plSsma confined' to a flux t^ which has substantial 


vari-atlon of magnetic field strength along its length. We shall concern 
ourselves with the phase in which the plasma cools primarily by thermal 
conduction to the chromosphere. This is a substantial simplification 
of the. real problem, which involves also radiation and mass motion. 
Nevertheless, we consider this study to be a useful first step in 
trying to understand the influence of magnetic-field geometry on the 
cooling of flare plasma. 

II. Model 

Since thermal conduction across magnetic-field lines is negligible 
in comparison with thermal conduction along field lines, we majr’ consider 
the cooling of a small flux tube independently of all other flux tubes. 
We adopt the model shown in Figure 1, showing a symmetrical flux tube 
which extends a height H above the top of the chromosphere which is 
the base of our model. The independent variable s measures distance 
along the loop from the top of the loop. Ihe variation of cross- 
sectional area is described by a function A(s), normalized to unity 
at the top of the loop. 

More specifically, we assume the field to be current-free above 
the chromosphere and to be identical with the field produced either by 
a horizontal line dipole (l) or by a hpriZdhtai point dipole (P) at 
a depth D below the chrpniosphere. It Is convenient to use the angle 
0 defined in Figure 1 in place of s. as the independent variable. The 
relationship between sand 0 is given by 


(2.1. P) 


I = I sin e (1+3 sln^e)^^^ + | sin 0 + (1+3 sin^e)^'^^] 

where 


R - H + D 


( 2 . 2 ) 


and we repeat equations for the line-dipole and point-dipole models, 
as necessary. The area function A(6), normalized to unity at 0 ^ o, 
is given by 


or 


A(6) = cos^0 

A(e) * cos^e (1 + 3 sin^e)"^^^ . 


(2.3.L) 


(2.3. P) 


For a plasma of temperature lo"^ K, the scale height is approximately 
the same as the solar radius. As a result, it Is a reasonable approximation 
to ignore gravitational effects. 

The key simplification in our model is to assume that each dependent 
variable may be expressed as a function of s multiplied by a function 
of t. In particular, we assume that the electron density n is expressible 
as 

n *'^‘'f (s) g(t) . (2.4) 

Zf, in addition, we assume that there Is no exchange of plasma between 
the flux tube and the chromosphere, it follows that g = constant so that 

n « n(s) . (2.5) 

It now ftjil'ows from the equatio^^^ of continuity that the plasma is static. 
Then, ignoring ^ravitatibnal effects, it follows from the equation of 


motion that the plaama must be isobaric, oo that 


P “ P(t) , 


( 2 . 6 ) 


where p is the total plasma pressure. For a fully ionized hydrogen 
plasma, p and n are related by 


p =< 2nkT 


(2.7) 


Since heating and radiation losses are being neglected, the heat 
equation is n imply 


|jC3nkI) . i |^(ak |I) 


( 2 . 8 ) 


where K is the themal conductivity. We adopt the Spitzer (1962) form. 


K - 


( 2 . 9 ) 


.-6 


where, with accuracy sufficient for our needs, we may take Of = 10 , 

On using equations (2.5), (2.6) and (2.7), we find that the heat 
equation may be rewritten in the following form. 


i n‘7/2 ^ . a 

P •’*■ Adards/ 


dt 


( 2 . 10 ) 


where we have introduced the symbol 


G(s) * I [2kn(s)] 


- 7/2 


( 2 . 11 ) 


It is to be noted that terms oh the left-hand side of equation (2.10) 
ate functions only of t.) whereas terms dn the right-hsiid side are f Unc tlona 

■ ■ ■ ■ .- V, V ■ ■ , , r . r. ■ 


only of s.^ Hence dach hide of the aqua tipri must be :eqa^^ cbrtstalqt 
which "we take to be ' * ^^bh this' choicei wb i f^nd n f it'Y'-'- 


^6 
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expressible as 


P(t) = p_(l + t/T)"^^^ , 


( 2 . 12 ) 


where p = p^ at t = 0. For the two cases considered, we find that 
G(0) = Gq - p"5/2 . 0 tan 9 


(c.l3*L) 


or 


where 


= °0 " IQ (2.X3.P) 


cos 6 


Yi(e) = (64 + 32 cos^e - IX cos*^0 - X5 cos^e) . (2.14) 

Since Yj^ (9) is a slowly* varying function of 9, confined to the range 
0.91 < Y < 1.09» we replace it by unity in subsequent formulas. 

We denote by 0j^ the value of 9 at the -”ba3e*' q £ the model identified 
with the chromosphere. This is related to the height of the loop and 
the depth of the dipole by 


or 


2 2 
H = R sin 0. , D = R cos 0, 


H ■= R(1 - cos^e) , D * R cos^e. 


( 2 .I 5 .L) 
(2.15. P) 


On noting that G(0jj) 0, we see from equations (2.13) that T is 
expressible as 


or 


where 


T » t* 9. tan 0. 
O b b 


t » i- tan^0j^ sec; 0," ,, 


^0 20 PQ 00 ^ 


( 2 .I 6 .L) 
( 2 .I 6 .P) 
(2.17) ^ 
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In this expression, » p(o) and » T(0,0). 

On using equations (2.?) through (2.I7), we obtain the following 
expression for the temperature: 

( cos G. 


2/7 


^ b 8 sin 9 

T(e,t) - T^Ct) 0 sin 8. cos 8 


(?.18.L) 


or 


TCe,t) = Tq(c) 1 - 


sin^ef^'^ 

2 I 

sin 0j^ cos 9 ) 


C 2 .I 8 .P) 


where T^ft), the temperature at the top of the loop, Is given by 

T(0,t) - T^Ct) = T^^(l + t/T)‘^^5 . 


(2.19) 


From equations (2.9) and (2.18), we find that the heat flux Is 
given by 

cos 0. 


F(6,t) - F„(t) 


O' ^ 0jj sin 0j^ 


6 + sin 9 cos 9 
2 cos^0 


(2.20.L) 


or 


where 


F(e,t) - F^(t) 


®b Sin 0 

0 '“" 2 — ^ Y2(0) 

^ sln^^G,. cos°e 


(2.20.P) 


Fg(t) » I Of (t) 


( 2 . 21 ) 


and 


^^(e) - (1 + 3 slh^0)^^^ (1 - sih®0 + |- sin^G - I s£n^9) . (2.2£») 

Ifc is fbund that YgC®) ® slowly yarding func|ti,on of 0, varying 

in the range ^94 < Yg 1 * 05 . Heofee in subseaueni qalculations we > 
replace y^ by unity^ . 


-- 




The heat flux F. (t) at the base of the tube Is found to be 


j 20jj ■*■ sln(2ej^) I * 


( 2 . 23 . 


Fb(t) « F^(t) 


sin $b cos ©b 


( 2 . 23 . 


llie total energy loss rate, referred to unit cross-sectional 
area at the top of the tube, is given by 


S = A(0b) F^Ct) 


(2.2lf) 


which leads to the express;ion 


i 2/i» 

cos 9. 


+ cot 0, 


( 2 . 25 . 


S( 0 b,t) - 2 FQ(t) cos^eb cot 6b (1 + 3 3in^6b)' 


(2.25. 


Ill . Cooling Rate 


in discussliVg r^)sults of our analysis of this model, it is convenient 

i • . * ' 

to Introduce the symbol T for the "compression factor" defined either as 
the area ratio 


( 34 ) 


or as the magnetic field ratio 



We aee from equation (2»3) that F is related to by 


r «■ sec^6. 

D 

or 

r o seAj^ (1 + 3 


(3.3.L) 


(3.3. P) 


We may now compare the cooling time t(F) of a loop structure with 
compression factor F to T„, the cooling time of a planar model with the 
same loop length Sj^. Since Tp is related to by 

Tp « Tq (Sj^/R)2 (3.h) 


we may calculate from equation (2.11) and (2.16) the dependence of 

T(F)Tp on F, as shown in Figure 2. We note that the cooling time is 

a sensitive function of F. The analysis by Rust and Bar (1973) of the 

loops produced by the flares of August, 1972* shew that, for those 

events, F took a range of values from 4 to 30. Such values clearly have 

an important effect in influencing heat loss. 

The initial; temperatures and densities of flare- produced soft 

x-ray-emltting regions are typically T = 10^*^ K and n cs 10^^ cm"^, 

with H ~ 10^*^ cm. On using equations (2.17) and (3.4), we find that 
8 

Tp ss 10 s. By comparison, t«v^ radiative cooling time may be estimated 
from 



a 

2 





where A ( 7 ) is the radiation loss function c a Icut a ted by Cox and Tucker 

(1969). We find,' fpr.^the pfoposed<'j>nnametefa> 10^*^ s. If 

the compression factor F is‘ large enough, the conduction cooling time 


may be increased by an order of magnitude to make it comparable with 
T„. In any case, radiation must bventually dominate since, as the 
plasma cools, increases as T decreases. Figure P. indicates that 
radiation will dominate at a higher temperature In a loop structure 
with large compression factor than it would in a loop structure with 
small compression factor. The fact that the influence of the magnetic- 
field geometry is to increase the cooling time, by comparison with a 
planar model, implies that a larger fraction of the initial energy 
will be radiated. Since it is unlikely that T will exceed 10^, it 
appears that the present model should be valid for the initial cooling 
phas<' of post-flare plasma loops. 

IV. Heat Flux and Evaporation 

Since the cooling rate decreases with increasing compression 
factor, it is clear that the total energy loss rate into the chromosphere 
must decrease as F increases. The ratio of the energy loss rate of a 
flare loop to that of a planar model is plotted against F in Figure 3, 
where we have used 

Fp^ Fo(t) . R/Sj, , (4.1) 

for the heat flux in the planar model. We see that the energy loss by 
thermal conduction may be decreased by an order of magnitude or more 
by the influence of magnetic-field georaet^, 

iSife ratio of the heat flux of a loop model to that of a planar 
model is shown in Figure 4 j from which we see that the heat , flux 
(per unit area) is greater for a Ippp structure than for a planar 
structure. :^e magtiebic field act4 as a futtnei increasing the heat 

10 


flux per unit area, but not rapidly enough to compensate for the 
decrease In area, so that the total energy loss by thermal conduction 
decreasrrs as the compression factor Increases. 

Even In the case of a planar model, one faces the difficulty that 
the high heat flux Into the chromosphere from a post- flare loop, of 
order 10^ erg cm s , cannot be maintained In the chromosphere 
because the thermal conductivity there is too low. Such fluxes are 
also top high to be converted into chromospheric radiation, it seems 
most likely that the excess heat flvnc into the chromosphere results 
in evaporation of chromospheric material. Since the heat fluxes are 
even higher when the effect of maghetlc-fleld geometry is taken into 
account, the likelihood of evaporation becomes greater. Evaporation 
will Increase the density in the loop, thus leadi.ng to a higher rate 
of energy loss by radiation. 

We recognize that a complete treatment of the evolution of fla're 
plasma must take actount of evaporation:/ We are developing the present 

. I ■ ■ 

model to Include this process. This more general study will help 
to delineate the range of applicability of the. present model, which 
ignores all plasma motion. Our present assessment is that evaporation 
will occur mainly at the /ihltial Stage in devaiopment of flare plasma, 


probably In re8pot;^$ii^ to bombardment of the"" chromosphere by a flux of . 
high“energy particles^ After the plasma Is evaporated^ there shbuld 
be ^ second stage described apprdxitaately by the pre model. After 
further cooling, radiation will beccN^i the dominant cooling S»chanlsi^ 


and thia ^preieht6:« tdiird .Stage idlih^^8 ^the~‘sub^^^ 




article being <prepared ^or publlcationv 








V* Dtncusslon 

The model developed la this article Is believed to be relevant 
to the early decay phase of poet-flare plasma. Although the model 
leads ~o a very simple form for the dependence of plasma temperature 
on time (equation (2.19)) • unlikely that this law can be 

compared directly with soft x-ray observations from post-flare plasma, 
since such plasma will typically occupy flux tube') covering a wide 
range of height and the decay constant V will vary from tube to tube. 
In order to test the model, observations of high spatial and temporal 
resolution are necessary so that both the temperature evolution and 
the geoDt''try of a flare loop can be determined. Alternatively, it 
would be possible to consider the behavior of -a model comprising a 
finite range of field lines and to integrate the total x-ray emission 
from this model. Such results might then be compared diroctly with 
observational data, 



The main aim of this article was to investigate the Influence 
upon flare Lodi Ing of realistic magnetic-field geometries. This is 
sunmari^ed by Figures 1 through b-,. In which the important parameter 
characterizing the geometry is the "compression factor" F, We see 


that magnetic-field geometry can hav^ a large effect on the flare 
behavior, indicating that it may be quite misleading to treat the 
evblutibn of flare plasma onMthe basis of a p line- par at lei model. 
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Figure Captions 

The geometry of our loop model, s is the distance along a 
field line sieasured from the top of the loop. A(a) is the 
croaa-sectlonal area of the flux tube. 

Ratio of cooling timed for a loop and for a planar (constant 
cross-section) model as a function of T, the ratio of the 
field strength at the base of the loop to that at the top. 

The solid line refers to a jfolnt dipole and the broken line 
to a line dipole. 

Ratio of energy loss rates for a loop and for a planar (constant 
cross-section) model as a. function of f, the ratio of the 
field strength at the base of the loop to thet at the top« 

The Solid line refers to a point dipole and the broken line 
to a line dipole'^ 

liatib of base hast fluxes for ,n loop and for a planar (cohstant 
cross-section) model as a function of T, the ratio of the 

<'» ■ '■ ' ' . ' - ''■■jv'''"'.''" 

■ ■ 1 * 1 . . 

field strength at the base of the loop to that at the top. 

The aolid line refer# to a ^olhf dipole and the broken line 

' ' ' , • ■ ' ■ a\ ■ ' ■ ■ 

■; to. a^-lln#' dipole*: ; 
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